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In an earlier paper we developed a Parton Bubble Model (PBM) for RHIC, high-energy heavy-ion 
collisions. PBM was based on a substructure of a ring of localized bubbles (gluonic hot spots) which 
initially contain 3-4 partons composed of almost entirely gluons. The bubble ring was perpendicular 
to the collider beam direction, centered on the beam, at midrapidity, and located on the expanding 
fireball surface of Au Au central collisions (0-10%) at y/sNN=200 GeV. The bubbles emitted cor- 
related particles at kinetic freezeout, leading to a lumpy fireball surface. For a selection of charged 
particles (0.8 GeV/c < pt < 4.0 GeV/c), the PBM reasonably quantitatively (within a few percent) 
explained high precision RHIC experimental correlation analyses in a manner which was consistent 
with the small observed HBT source size in this transverse momentum range. We demonstrated 
that surface emission from a distributed set of surface sources (as in the PBM) was necessary to 
obtain this consistency. In this paper we give a review of the above comparison to central Au Au 
collisions. The bubble formation can be associated with gluonic objects predicted by a Glasma Flux 
Tube Model (GFTM) that formed longitudinal fiux tubes in the transverse plane of two colliding 
sheets of Color Glass Condensate (CGC), which pass through one another. These sheets create 
boost invariant fiux tubes of longitudinal color electric and magnetic fields. A blast wave gives the 
tubes near the surface transverse fiow in the same way it gave transverse fiow to the bubbles in 
the PBM. In this paper we also consider the equivalent characteristics of the PBM and GFTM and 
connect the two models. In the GFTM the longitudinal color electric and magnetic fields have a 
non-zero topological charge density FF. These fields cause a local strong CP violation which effects 
charged particle production coming from quarks and anti-quarks created in the tube or bubble. 

PACS numbers: 25.75.Nq, ll.30.Er, 25.75.Gz, 12.38.Mh 



I. PARTON BUBBLE MODEL 
DEVELOPMENT 



Our interest in and the eventual development of 
the PBM goes back many years to the early nineteen 
eighties. Van Hove's work[l| in the early eighties 
considered the bubbles as small droplets of quark- 
gluon plasma (QGP) assumed to be produced in 
ultra-relativistic collisions of hadrons and/or nu- 
clei. His work recognized that experimental detec- 
tors only see what is emitted in the final state at 
kinetic freezeout. Therefore final state predictions 
of his string model work on plasma bubbles or any 
other theory required specific convincing predictions 
that are observable in the final state which can be 
measured by experiment. His work predicted that 
the final state rapidity distribution dn/dy of hadrons 
would exhibit isolated maxima of width Ay ^ 1 (sin- 
gle bubble) or rapidity bumps of a few units (due to 
those cases producing a few bubbles). These rapid- 
ity regions would also contain traditional signals of 
QGP formation. Even this early work by Van Hove 
needed his bubbles near the surface of the final state 
fireball in order to experimentally measure the par- 
ticles emitted in the final state at kinetic freezeout. 
We and other searched for these Van Hove bubbles, 
but no one ever found any significant evidence for 



them. 

Our next paper Q was our final attempt to make a 
theoretical treatment of the single bubble case (sim- 
ilar to the Van Hove case) which could be experi- 
mentally verified. It is possible that with enough 
statistics one in principle could find single to a few 
bubble events. We developed a number of event gen- 
erators which could possibly provide evidence for 
striking signals resulting from these bubbles. One 
should note that in that paper we had included in 
Sec. 11 mathematical expressions describing how 
charged pions are effected by strong color electric 
and color magnetic fields which are present and par- 
allel in a CP-odd bubble of metastable vacuum. This 
section was motivated by the work of Kharzeev and 
Pisarski 3]. For example eq. 3-4 of the section show 
the boosts of the 7r+ and the tt~ we estimated. Thus 
in 2000 we were already interested in and were in- 
vestigating and publishing predictions for this phe- 
nomenon. 

When we considered the RHIC quantum interfer- 
ence data in 20000 it became clear that the fire- 
ball surface was rapidly moving outward. This im- 
plied a very large phase space region covered by the 
fireball, and also implied that it would be very un- 
likely for there to be only one isolated bubble (small 
source size) with a large amount of energy sitting 
on the surface. It seemed more likely that there 
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would be many bubbles or gluonic hot spots around 
the expanding surface. This led to a paper 5] which 
was our earlier version of Ref.(6||. In that paper 
we concluded that the behavior of the Hanbury- 
Brown and Twiss (HBT) measurements [1, 0] should 
be interpreted as evidence of a substructure of bub- 
bles located on the surface of the final state fireball 
in the central rapidity region at kinetic freezeout. 
The HBT radii were decreasing almost linearly with 
transverse momentum and implying to us a source 
size of ~2 fm radii bubbles were on the surface and 
could be selected for if we considered transverse mo- 
menta above 0.8 GeV/c. These momenta would al- 
low sufficient resolution to resolve individual bubbles 
of ^2 fm radii. We further concluded that these 
HBT quantum interference observations were likely 
due to phase space focusing of the bubbles pushed 
by the expanding fireball. Thus the HBT measure- 
ments of source size extrapolated to transverse mo- 
menta above 0.8 GeV/c were images of the bubbles. 
The HBT correlation has the property of focusing 
these images on top of each other for a ring of bub- 
bles transverse to and centered on the beam forming 
an average HBT radius. 



Thus our model was a ring of bubbles sitting on 
the freezeout surface with average size ~2 fm ra- 
dius perpendicular to the beam at mid-rapidity. The 
phase space focusing would also lead to angular cor- 
relations between particles emitted by the bubbles 
and be observable. Fig. 1 shows the geometry of the 
bubbles. For the background particles that account 
for particles in addition to the bubble particles we 
used unquenched HIJING[1] with its jets removed. 
We assumed that most of the jets were eliminated 
from the central events because of the strong jet 
quenching observed at RHIC^Q,]- We investigated the 
feasibility of using charged-particle-pair correlations 
as a function of angles which should be observable 
due to the phase space focusing of the particles com- 
ing from individual bubbles. At that point we fully 
expected that these correlations would be observable 
in the STAR detector 10] at RHIC if one analyzed 
central Au Au collisions. Correlation analyses are 
powerful tools in detecting substructures. Histor- 
ically substructures have played an important role 
in advancing our understanding of strong (non per- 
turbative) interactions. Our earlier paper explained 
the general characteristics of the angular correlation 
data, and was also consistent with HBT measure- 
ments in a qualitative manner. This motivated us to 
develop a reasonably quantitative model, the parton 
bubble modeljsl which is discussed in the following. 



II. PARTON BUBBLE MODEL [6] 

In this publication 'g'] we developed a QCD in- 
spired parton bubble model (PBM) for central (im- 
pact parameter near zero) high energy heavy ion col- 
lisions at RHIC. The PBM is based on a substruc- 
ture consisting of a single ring of a dozen adjoining 
2-fm-radius bubbles (gluonic hot spots) transverse to 
the collider beam direction, centered on the beam, 
and located at or near mid-rapidity. The at or near 
mid-rapidity refers to the boost invariant region that 
exist in the RHIC heavy ion collisions which spans 
2-4 units of rapidity. Because of this spread in rj, 
we will use the relative spread of particles in 77 [Ar]) 
in order to capture the angular correlations along 
the beam axis. The ring resides on the fireball blast 
wave surface (see Fig. 1). We assumed these bub- 
bles are likely the final state result of quark-gluon- 
plasma (QGP) formation since the energy densities 
produced experimentally are greater than those es- 
timated as necessary for formation of a quark-gluon- 
plasma. Thus this is the geometry for the final state 
kinetic freezeout of the QGP bubbles on the sur- 
face of the expanding fireball treated in a blast wave 
model. 

The average behavior of emitted final state par- 
ticles coming from the surface bubbles at kinetic 
freezeout is given by the ring of twelve bubbles 
formed by energy density fluctuations near the sur- 
face of the expanding fireball of the blast wave. One 
should note that the blast wave surface is moving at 
its maximum velocity at freezeout (3c/4). For cen- 
tral events each of the twelve bubbles have 3-4 par- 
tons per bubble each at a fixed (j) for a given bubble. 
This number of partons was determined from cor- 
relation data analyzed by the STAR experiment @. 
The transverse momentum (pt) distribution of the 
charged particles is similar to pQCD but has a sup- 
pression at high pt like the data. 

The bubble ring radius of our model was estimated 
by blast wave, HBT and other general considerations 
to be approximately 8 fm. The bubbles emit corre- 
lated charged particles at final-state kinetic freeze- 
out where we select a pt range (0.8 GeV/c < pt < 
4.0 GeV/c) in order to increase signal to background. 
The 0.8 GeV/c pt cut increases the resolution to al- 
low resolving individual bubbles which have a radius 
of ~2 fm. This space momentum correlation of the 
blast wave provides us with a strong angular corre- 
lation signal. PYTHIA fragmentation functions [ll| 
were used for the bubbles fragmentation that gener- 
ate the final state particles emitted from the bubbles. 
A single parton using PYTHIA forms a jet with the 
parton having a fixed r] and (f) (see Fig. 2). The 3-4 
partons in the bubble which shower using PYTHIA 
all have a different 77 value but all have the same </> 
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Plane section of bubble geometry perpendicular to the beam at eta = 




A perspective view of the bubble geometry. 




Beam axis 



FIG. 1: The bubble geometry is an 8 fm radius ring perpendicular to and centered on the beam axis. It is composed 
of twelve adjacent 2 fm radius spherical bubbles elongated along the beam direction by the Landau longitudinal 
expansion. The upper left figure is a projection on a plane section perpendicular to the beam a:xis. The lower 
left figure is a projection of the bubble geometry on a plane containing the beam axis. The lower right figure is a 
perspective view of the bubble geometry. Due to rotational invariance about the beam axis the only direction that is 
meaningful to define is the beam axis shown in the lower 2 figures. 



(see Fig. 3). The PBM explained the high precision 
Au Au central (0-10%) collisions at ^/snn = 200 
GeV[12,] (the highest RHIC energy). 



A. The Correlation Function 

We utilize a two particle correlation function in 
the two dimensional (2-D) space of A<j> vs A77. The 
azimuthal angle of a particle is defined by the angle 
of the particle with respect to the vertical axis which 
is perpendicular to the beam axis and is measured 
in a clock-wise direction about the beam. A0 is the 
difference, - 4>2, of the </> angle of a pair of particles 
(1 and 2). The pseudo-rapidity 77 of a particle is 
measured along one of the beam directions. A77 is 
the difference, 771 - 772, of the rj values of a pair of 
particles (1 and 2). 

The two dimensional (2-D) total correlation func- 
tion is defined as: 



C(A^, A77) = S{Acl>, A77)/M(A0, A77). (1) 



Where S(A(/), A?]) is the number of pairs at the 
corresponding values of A0, Arj coming from the 
same event, after we have summed over all the 
events. M(A0, A77) is the number of pairs at the cor- 
responding values of A^, Ar/ coming from the mixed 
events, after we have summed over all our created 
mixed events. A mixed event pair has each of the 
two particles chosen from a different event. We make 
on the order of ten times the number of mixed events 
as real events. We rescale the number of pairs in the 
mixed events to be equal to the number of pairs in 
the real events. This procedure implies a binning 
in order to deal with finite statistics. To enhance 
comparison we use the same binning in our simula- 
tions as is used by the STAR high precision experi- 
mental analyses 12|. The division by M{A(j>, Arj) for 
experimental data essentially removes or drastically 
reduces acceptance and instrumental effects. If the 
mixed pair distribution was the same as the real pair 
distribution C{A(j), A77) would be one for all values of 
the binned A0, Aij. In the correlations used in this 
paper we select particles independent of its charge. 
The correlation of this type is called a Charge In- 
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FIG. 2: A jet parton shower. 



dependent (CI) Correlation. This difference correla- 
tion function has the defined property that it only 
depends on the differences of the azimuthal angle 
(A(/)) and the beam angle (A77) for the two particle 
pair. Thus the two dimensional difference correla- 
tion distribution for each bubble which is part of 
C{Aip, A77) is similar for each of our 12 bubbles and 
will image on top of each other. The PBM fit to 
the angular correlation data was reasonably quanti- 
tative to within a few percent (see Fig. 4 and Fig. 
5)). 

The correlation functions we employed (like the 
HBT correlation functions) have the property that 
for the difference in angles (difference of momentum 
for HBT) these correlations will image all 12 bub- 
bles on top of each other. This leads to average 
observed angles of approximately 30° in A0, 70° 
in A77, originating from a source size of about ^2 
fm radius which is consistent with the HBT correla- 
tion. Thus the PBM generates A(f) vs A?/ charged- 
particle-pair correlations for charged particles with 



Pt in the range 0.8 Gev/c to 4.0 GeV/c as displayed 
in Fig. 4 and Fig. 5. We show the CI correla- 
tion in Fig. 5 a reasonably quantitative successful 
comparison with data. It should be noted that the 
spread in A77 of the bubble is determined by the 
longitudinal momentum distribution of the 3-4 par- 
tons that make up the bubble. This distribution 
which is quassian in nature was adjusted to the ob- 
served Ar] width of correlation data analyzed by the 
STAR experiment fi]. Furthermore the model results 
were consistent with the Hanbury-Brown and Twiss 
(HBT) observations 7] that the observed source radii 
determined by quantum statistical interference were 
reducing by a considerable factor with increasing 
transverse momentum (pt)- The HBT radii are in- 
terpreted to reduce from ^6 fm at pt ~0.2 GeV/c 
to ^2 fm at Pt ^1 GeV/c for our pt range. The gen- 
erally accepted explanation for this behavior is that 
as Pt increases radial flow increasingly focused the 
viewed region of the final state into smaller volumes. 
If just one small HBT size bubble were emitting all 
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FIG. 3: Each bubble contains 3-4 partons as shown. 



the correlated particles, this phenomena would lead 
to large spikes of particles emitted at one limited 
(j) angular region in individual events. This is defi- 
nitely not observed in the Au Au or other collision 
data at RHIC. Therefore, a distributed ring of small 
sources around the beam as assumed in the PBM is 
necessary to explain both the HBT results [7*1 and the 
correlation data[12]. The particles emitted from the 
same bubble are virtually uncorrelated to particles 
emitted from any other bubble except for momen- 
tum conservation requirements. An away side peak 
in the total correlation is built up from momentum 
conservation between the bubbles. 



III. THE PBM FOR OTHER 
CENTRALITIES(PBME[l3l) 

The PBM was also recently extended to 
PBME[13| which is identical to the PBM for cen- 
tral collisions (0-5%). For centralities running from 



30-80% jet quenching is not strong enough to make 
jets negligible therefore, a jet component was added 
which was based on HIJING calculations. This jet 
component accounts for more of the correlation as 
one moves toward peripheral bins, and explains all of 
it for the most peripheral collisions. The PBME ex- 
plained in a reasonably quantitative manner (within 
a few percent) the behavior of the recent quantita- 
tive experimental analysis of charge pair correlations 
as a function of centrality fl^ . This further strength- 
ened the substantial evidence for bubble substruc- 
ture. The agreement of the PBM and the PBME 
surface emission models with experimental analyses 
strongly implied that at kinetic freezeout the fire- 
ball was dense and opaque in the central region and 
most centralities (except the peripheral region) in 
the intermediate transverse momentum region (0.8 
< Pt < 4.0 GeV/c). Thus we conclude that the ob- 
served correlated particles are both formed and emit- 
ted from or near the surface of the fireball. In the 
peripheral region the path to the surface is always 



FIG. 4: The CI (sum of aU charged-particle-pairs) correlation for the 0-10% centrality bin with a charged particle pt 
selection (0.8 < pt < 4.0 GeV/c), generated by the PBM. It is plotted as a two dimensional A(j> vs. Arj perspective 
plot. 



small. 

The CI correlation used in Fig. 4 and the PBM [6] 
is also defined as the sum of two charged-particle- 
pair correlations. They are the unlike-sign charge 
pairs (US) and the like-sign charge pairs (LS). This 
Charge Independent correlation (CI) is defined as 
the correlation made up of charged-particle-pairs in- 
dependent of what the sign is, which is the average 
of the US plus the LS correlations. Thus CI = (US 
-I- LS)/2. Therefore the CI is the total charge pair 
correlation observed in the experimental detection 
system within its acceptance. 

The CI was generated for ^/sNN = 200 GeV 
central (0-5%) Au An collisions in the pt range 
0.8 GeV/c to 4.0 GeV/c at RHIC. When the CI 
was corrmared to the corresponding experimental 
analyses (l2i a reasonably quantitative agreement 
within a few percent of the observed CI was at- 
tained. See Ref.0 Sec. 4.2, RehfTsj Sec. V, and 
Ref.[lJ] Sec. VI B^. This analysis and our previous 
workQ was done without using a jet trigger. In 
fact jet production in the central region was negligi- 



^ The STAR collaboration preferred to use the conventional 
definition of CI = US -I- LS which is larger by a factor of 
2 than the CI definition used in this present paper that is 
directly physically meaningful. 



ble due to strong jet quenching [ol. [Tsl. [l6j . As one can 
see in Fig. 4 the central CI correlation of Acj) on the 
near side (At/) < 90°), is sharp and approximately 
jetlike at all Ary. In contrast the Arj dependence is 
relatively flat. In the central production the average 
observed angles are approximately 30° in A(j) and 
70° in A77. It is of interest to note that the sharp 
jetlike collimation in A(/) persists at all centralities 
(0-80%). The elongation in Arj persists in the 0-30% 
centrality range and then decreases with decreasing 
centrality becoming jetlike in the peripheral bins. 
These characteristics are in reasonable quantitative 
agreement with the fits of the PBM and the PBME 
to experimental data analyses (gI [T2| - [i3 1 . 



The difference of the US and LS correlations is 
defined as the Charge Dependent (CD) correlation 
(CD = US - LS). The 2-D experimental CD cor- 
relation for centralities (0-80%) has a jetlike shape 
which is consistent with PYTHIA jets (vacuum frag- 
mentation) for all centralities. This clearly implies 
that both particle hadronization and emission oc- 
cur from the fireball surface region as explained in 
Ref.[ll Sec. II. 
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FIG. 5: In each of the five labeled A77 bins we show the A(/) total correlation for the CI as a function of A^. The 
STAR Au + Au central trigger analysis results from the formulae of Ref. [l^ are shown as a solid line. The parton 
bubble model predictions are shown by the circular points(o) which are large enough to include the statistical errors 
from a 2 million event sample. The vertical correlation scale is not offset and is correct for the largest A?7 bin, 1.2 
< A?7 < 1.5, and is the lowest bin on the figure. As one proceeds upwards to the next bin A77 the correlation is offset 
by +0.01. This is added to the correlation of each subsequent A»7 bin. The smallest A?7 and top of Figure 12 has a 
+0.04 offset. A solid straight horizontal line shows the offset for each Ar; bin. Each solid straight horizontal line is 
at 1.0 in correlation strength. The 0-10% centrality in HIJING corresponds to impact parameter (b) range 0.0 to 4.0 
fm. The agreement is very good. 



IV. GLASMA FLUX TUBE MODEL 

A glasma flux tube model (GFTM)[17] that had 
been developed considers that the wavefunctions of 
the incoming projectiles, form sheets of color glass 
condensates (CGC)[3 that at high energies colhde, 
interact, and evolve into high intensity color elec- 
tric and magnetic fields. This collection of primor- 
dial fields is the Glasma [H, [13], and initially it is 
composed of only rapidity independent longitudinal 
color electric and magnetic fields. An essential fea- 
ture of the Glasma is that the fields are localized 
in the transverse space of the collision zone with a 



size of 1 IQs- Qs is the saturation momentum of par- 
tons in the nuclear wavefunction. These longitudinal 
color electric and magnetic fields generate topologi- 
cal Chern-Simons charge [2l| which becomes a source 
for particle production. 

The transverse space is filled with flux tubes of 
large longitudinal extent but small transverse size 
^Q^^- Particle production from a flux tube is a 
Poisson process, since the flux tube is a coherent 
state. As the partons emitted from these flux tubes 
locally equilibrate, transverse flow builds due to the 
radial flow of the blast wave[22j. The flux tubes that 
are near the surface of the fireball get the largest 



8 



radial flow and are emitted from the surface. As in 
the parton bubble model these partons shower and 
the higher pt particles escape the surface and do not 
interact. These flux tubes are in the present paper 
considered strongly connected to the bubbles of the 
PBM. The method used to connect the PBM and 
the GFTM is described and discussed in the next 
section V. 

Qs is around 1 GeV/c thus the transverse size of 
the flux tube is about 1/4 fm. The flux tubes near 
the surface are initially at a radius ~5 fm. The </> an- 
gle wedge of the flux tube is ~l/20 radians or ~3°. 
Thus the flux tube initially has a narrow range in 
(f). The large width in the Ar] correlation which in 
the PBM depended on the large spread in A77 of the 
bubble partons results from the independent longitu- 
dinal color electric and magnetic fields that created 
the Glasma flux tubes. How much of these longitudi- 
nal color electric and magnetic fields are still present 
in the surface flux tubes when they have been pushed 
by the blast wave will be a speculation of this paper 
for measuring strong CP violation? 

It has been noted that significant features of the 
PBM that generated final state correlations which fit 
the experimentally observed correlation datafs', 
[ll] are similar to those predicted by the GFTM 17]. 
Therefore in the immediately following section we 
assume a direct connection of the PBM and the 
GTM, give reasons to justify it, and then discuss 
its consequences, predictions and successes. 



V. THE CONNECTION OF THE PBM AND 
THE GFTM 

The successes of the PBM have strongly implied 
that the final state surface region bubbles of the 
PBM represent a significant substructure. In this 
subsection we show that the characteristics of our 
PBM originally developed to fit the precision STAR 
Au Au correlation data in a manner consistent with 
the HBT data; implies that the bubble substructure 
we originally used to fit these previous data is closely 
related to the GFTM. 

A natural way to connect the PBM bubbles to the 
GFTM fiux tubes is to assume that the final state 
at kinetic freezeout of a flux tube is a PBM final 
state bubble. Thus the initial transverse size of a 
flux tube ^1/4: fm has expanded to the size of ^2 
fm at kinetic freezeout. With this assumption we 
find consistency with the theoretical expectations of 
the GFTM. We can generate and explain the trig- 
gered ridge phenomenon and data (see Sec. VI), thus 
implying the ridge is connected to the bubble sub- 
structure. We can predict and obtain very strong ev- 
idence for the color electric field of the glasma from 



comparing our multi-particle charged particle cor- 
relation predictions, and existing experimental cor- 
relation publications (see Sec. VII). We have also 
predicted correlations which can be used to search 
for evidence for the glasma color magnetic field (Sec. 
VII). 

Of course an obvious question that arises is that 
since a fiux tube is an isolated system does a PBM 
bubble, that we assume is the final state of a 
flux tube at kinetic freezeout, also act as an iso- 
lated system when emitting the final state parti- 
cles? The near side correlations signals since the 
original PBmQ have always come virtually entirely 
from particles emitted from the same bubble. Thus 
each bubble has always acted as an isolated system 
similar to the behavior of a fiux tube. 



VI. THE RIDGE IS FORMED BY THE 
BUBBLES WHEN A JET TRIGGER IS 
ADDED TO THE PBM 

In heavy ion collisions at RHIC there has been 
observed a phenomenon called the ridge which has 
many different explanations [TtI [23l - l30| . The ridge is 
a long range charged particle correlation in A77 (very 
flat), while the A4> correlation is approximately jet- 
like (a narrow Gaussian). There also appears with 
the ridge a jet-like charged-particle-pair correlation 
which is symmetric in Ary and A0 such that the peak 
on the jet-like correlation is at Arj = and A(/) = 0. 
The A(j) correlation of the jet and the ridge are ap- 
proximately the same and smoothly blend into each 
other. The ridge correlation is generated when one 
triggers on an intermediate pt range charged particle 
and then forms pairs between that trigger particle 
and each of all other intermediate charged particles 
with a smaller pt down to some lower limit. The first 
case we will study in this paper is a trigger charged 
particle between 3.0 to 4.0 GeV/c correlated with all 
other charged particles which have a pt between 1 . 1 
GeV/c to 3.0 GeV/c. 

In this paper we will investigate whether the PBM 
can account for the ridge once we add a jet trigger 
to our PBM generator j^. However this trigger will 
also select jets which previously could be neglected 
because there was such strong quenching d, [ll, [l^ 
of jets in central collisions. A jet trigger had not 
been used in the PBM comparison to all previous 
data. We use HIJING^Sl merely to determine the 
expected number of jets to add for our added jet 
trigger. These jet particles were added to our PBM 
generator. Thus our PBM generator now had HI- 
JING generated background particles, bubbles of the 
PBM and added jet particles from HIJING. We have 
already shown that our final state particles come 
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from hadrons at or near the fireball surface. We re- 
duce the number of jets by 80% which corresponds 
to the estimate that only the parton interactions on 
or near the surface are not quenched away, and thus 
at kinetic freezeout form and emit hadrons which 
enter the detector. This 80% reduction is consistent 
with single 7r° suppression observed in Ref . [l6| . We 
find for the reduced HIJING jets that 4% of the Au 
Au central events (0-5%) centrality at ^/SNN — 200 
have a charged particle with a pt between 3.0 and 4.0 
GeV/c with at least one other charged particle with 
its Pt greater than 1.1 GeV/c coming from the same 
jet. The addition of the jets to the PBM generator 
provides the appropriate particles which are picked 
up by the trigger in order to form a narrow Arj cor- 
relation signal at which is also a narrow signal in 
Aip at (Fig. 14). This narrow jet signal is present 
in the data and is what remains of jets after 80% are 
quenched away. 

We then form two-charged-particle correlations 
between one-charged-particle with a pt between 3.0 
to 4.0 GeV/c and another charged particle whose 
Pt is greater than 1.1 GeV/c. The results of these 
correlations are shown in Fig. 6. Fig. 6 is the CI 
correlation for the 0-5% centrality bin with the just 
above stated pt selections on charge pairs. Since we 
know in our Monte Carlo which particles are emit- 
ted from bubbles and thus form the ridge, we can 
predict the shape of the ridge for the above pt cut 
by plotting only the correlation formed from pairs of 
particles that are emitted by the same bubble (see 
Fig. 7). The charge pair correlations are virtually 
all emitted from the same bubble. Those charge 
pair correlations formed from particles originating 
from different bubbles are small contributors which 
mainly have some effect on the away side correla- 
tion. Thus Fig. 7 is the ridge signal which is the 
piece of the CI correlation for the 0-5% centrality 
of Fig. 6, after removing all other pairs except the 
pairs emitted from the same bubble. 



A. Parton pt correlated vs. random 

A very important aspect of the GFTM is the boost 
that the flux tubes get from the radial flow of the 
blast wave. This boost is the same all along the flux 
tube and only depends on how far away from the 
center axis of the blast wave the flux tube is. In 
the PBM the partons in the bubble received a boost 
in Pt from the radial flow field in the blast wave; 
that after final state fragmentation of the bubbles 
gave the generated particles a consistent pt spectrum 
with the data. Since the boost from the blast wave 
depended on the position of the bubble in the ra- 
dial flow field, there should be a correlation between 



partons pt within a bubble. If one redistributed the 
partons with their pt boosts uniformly among the 
bubbles, the over all results for the generated corre- 
lations without a trigger would be unchanged. One 
should note none of our prior work (i.e. PBM 6] or 
PBME[13]) contained a trigger. We have only added 
a trigger here to treat the "triggered ridge" which 
obviously requires it. In our treatment of the GFTM 
(Sec. VII) we have removed the trigger. 

Once we require a trigger demanding higher pt 
particles, we start picking out bubbles which have 
more radial flow (harder particles). Thus correlated 
particles which pass our pt selection and trigger come 
almost entirely from the harder bubbles with more 
radial flow, while the softer bubbles which were sub- 
ject to less radial fiow mainly generate low enough 
Pt particles that become background particles or are 
lost to the correlation analysis due to pt selection. 
Figure 6 shows the result of our trigger and pt se- 
lection. If one redistributed the partons with their 
Pt radial flow boost among all the bubbles, then all 
bubbles become equal and there is no longer soft 
and hard bubbles. With this change we can generate 
bubble events. We find that both the non-triggered 
correlation and the correlation are virtually un- 
changed in the region less than about 120° (see Fig. 
8). It seems that we move from the situation where 
we have a few bubbles with a lot of correlated parti- 
cles above our pt cut to a lot of bubbles that have few 
correlated particles above our cuts. However there is 
a difference in the triggered correlation's away side 
or A(j) near the 180° peak. 

In Ref.Q] we discuss this away side effect and at- 
tributed this peak to momentum conservation. This 
is however not the whole story. The away side peak 
is generated by the geometry of the bubble ring 
which requires that for every bubble there is another 
bubble nearly 180° away that is emitting particles. 
These symmetry requirements on the geometry con- 
tribute more to the away side peak especially for 
the case without a trigger or redistributed random 
partons. 

This geometry effect is what causes the away side 
correlation of elliptic fiow in a peripheral heavy ion 
collision. In such a collision there is a higher en- 
ergy density aligned with the reaction plane. On 
one side of the reaction plane there are more par- 
ticles produced because of increased energy density, 
and because of the geometry of the situation. On 
the other side of the reaction plane there are also 
more particles produced for the same reasons. Thus 
there is a correlation between particles that are near 
each other on one side of the reaction plane and a 
correlation between particles that are on both sides 
of the reaction plane. Momentum conservation re- 
quires particle communication by bouncing into each 
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Triggered particles 




FIG. 6: The CI correlation for the 0-5% centrahty bin that results from requiring one trigger particle pt above 3 
GeV/c and another particle pt above 1.1 GeV/c. It is plotted as a two dimensional A(jf> vs. A77 perspective plot. 



other, while geometry has no such communication. 
The right hand does not need to know what the left 
hand is doing. However, if physics symmetry makes 
them do the same thing then they can show a cor- 
relation. 



B. The away side peak 

The away side peak depends on the fragmenta- 
tion of the away side bubble. For the triggered case 
where there is a hard bubble with a strong corre- 
lation of parton pt inside the bubble (GFTM like), 
we will have more correlated particles adding to the 
signal on the near side. On the other hand the away 
side bubble will be on average a softer bubble which 
will have a lot less particles passing the pt cuts and 
thus have a smaller signal. If we consider the case 
of a random parton pt bubble the particles passing 
the Pt cuts should be very similar on the near and 
the away side. In Fig. 8 we compare the two cor- 
relations generated by the correlated and the ran- 
dom cases. We compare the A0 correlation for the 
two cases (solid is correlated and open is random) for 
each of the five Ayy bins which cover the entire Ar? 
range (0.0 to 1.5). The vertical correlation scale is 
not offset and is correct for the largest Ar/ bin, 1.2 
< A?7 < 1.5, which is the lowest bin on the figure. 
As one proceeds upwards to the next A77 bin the 



correlation is offset by -1-0.05. This is added to the 
correlation of each subsequent A77 bin. The smallest 
Ary bin on top of Fig. 8 has a -1-0.2 offset. A solid 
straight horizontal line shows the offset for each Ary 
bin. Each solid straight horizontal line is at 1.0 in 
correlation strength. We see that the away side cor- 
relation in the A0 region greater than about 120° is 
larger for the random case. 



C. Predicted transverse momentum 
dependence of correlation 

Particles which come from bubbles that satisfy the 
trigger (3.0 < pt < 4.0 GeV/c) are harder and have a 
different pt distribution compared to the case with- 
out trigger requirements. In a given triggered event 
if one could select particles that come from the bub- 
ble or ridge, one would find that the pt spectrum is 
harder than the average pt spectrum. In Fig. 9 we 
show how one can by trigger choose particles that 
are rich with ridge particles. We consider a typi- 
cal trigger particle and the associated regions in the 
available experimental (j) and rj ranges. The ridge 
particles we consider are separated from the trigger 
particle by 0.7 in rj. This is done to eliminate par- 
ticles coming from jet production that could also be 
associated with the trigger particle. The bulk of the 
ridge particles lie within 20° of the (p of the trigger 
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Ridge signal 




FIG. 7: The ridge signal is the piece of the CI correlation for the 0-5% centrality of Fig. 6 after removing all other 
particle pairs except the pairs that come from the same bubble. It is plotted as a two dimensional l\ij> vs. Ar; 
perspective plot. As explained in the text charge pair correlations formed by particles emitted from different bubbles 
are small contributors which mainly have some effect on the away side correlation for greater than about 120° . 



particle. The width of the spread is A0 = 40°. 
Therefore on a trigger by trigger basis we accumu- 
late particles from the ridge region (see Fig. 9) and 
form a pt spectrum of charged particles. We then 
form a ratio of the pt spectrum for the ridge region 
particles divided by the pt spectrum for all parti- 
cles in central Au Au events. This ratio is highly 
dependent on our ridge cut area. We can virtually 
remove this dependency by normalizing the two pt 
spectra to have the same counts for the 2.0 GeV/c 
bin. Figure 10 is this ratio which can be considered 
a correlation of particles as a function of pt in the 
ridge region compared to particles from the event in 
general. In Fig. 10 we form this ratio or correlation 
for the bubbles which have a correlated pt among 
the partons (GFTM like). We also show the correla- 
tion for bubbles which have a random pt distribution 
among partons in the bubble. For the random case 
the Pt spectrum does not differ much from the aver- 
age Pt spectrum while there is a big difference for the 
flux tube like case. This correlation could be easily 
measured in the RHIC data. 



D. Comparison to data 



Triggered angular correlation data showing the 
ridge was presented at Quark Matter 2006 [3l|. Fig- 
ure 11 shows the experimental A0 vs. A77 CI cor- 
relation for the 0-10% central Au Au collisions at 
y/SNN — 200 GeV; requiring one trigger particle pt 
between 3 and 4 GeV/c and an associated particle 
Pt above 2.0 GeV/c. The yield is corrected for the 
finite A77 pair acceptance. For the PBM generator, 
we then form a two-charged-particle correlation be- 
tween one charged particle with a pt between 3.0 to 
4.0 GeV/c and another charged particle whose pt 
is greater than 2.0 GeV/c. The results of this cor- 
relation is shown in Fig. 12. Figure 11 shows the 
corrected pair yield determined in the central data 
whereas Fig. 12 shows the correlation function gen- 
erated by the PBM which does not depend on the 
number of events analyzed. We can compare the 
two figures, if we realize that the away side ridge has 
around 420,000 pairs in Fig. 11 while in Fig. 12 the 
away side ridge has a correlation of around 0.995. 
If we multiply the correlation scale of Fig. 12 by 
422,111 in order to achieve the number of pairs seen 
in Fig. 11, the away side ridge would be at 420,000 
and the peak would be at 465,000. This would make 
a good agreement between the two figures. 
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FIG. 8: The A<j!> CI correlation for the 0-5% centrality bin requiring one trigger particle pt above 3 GeV/c and another 
particle pt above 1.1 GeV/c. We compare the two correlations generated by the correlated and the random cases. 
We compare the A(j> correlation for the two cases (solid is correlated and open is random) for each of the five Ar; bins 
which cover the entire A-q range (0.0 to 1.5). The vertical correlation scale is not offset and is correct for the largest 
Arj bin, 1.2 < A77 < 1.5, which is the lowest bin on the figure. As one proceeds upwards to the next Arj bin the 
correlation is offset by +0.05. This is added to the correlation of each subsequent Arj bin. The smallest Arj bin on 
top of Fig. 8 has a +0.2 offset. A solid straight horizontal line shows the offset for each Arj bin. Each solid straight 
horizontal line is at 1.0 in correlation strength. We see that the away side correlation for A<j) greater than about 120° 
is larger for the random case. 



The correlation formed by the ridge particles is 
generated almost entirely by particles emitted by 
the same bubble. We have shown in all our pub- 
lications that the same side correlation signals are 
almost entirely formed by particles coming from the 
same bubble. Thus we can predict the shape and 
the yield of the ridge for the above pt trigger selec- 
tion and lower cut, by plotting only the correlation 
coming from pairs of particles that are emitted by 
the same bubble (see Fig. 13). 

In Ref . it was assumed that the ridge yield was 
flat across the acceptance while in Fig. 13 we see 
that this is not the case. Therefore our ridge yield is 



approximately 35% larger than estimated in Ref. [3l| . 
Finally we can plot the jet yield that we had put into 
our Monte Carlo. We used HIJING^ to determine 
the number of expected jets, and then reduced the 
number of jets by 80%. This assumes that only the 
parton interactions on or near the surface that form 
hadrons at kinetic freezeout are not quenched away 
and thus enter the detector. This 80% reduction 
is consistent with single 7r° suppression observed in 
Ref.[16i]. This jet yield is plotted in Fig. 14 where 
we subtracted contributions from the bubbles and 
the background particles from Fig. 11. 
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FIG. 9: In this figure wc display the d> and r; ranges considered for producing ridge particles in the Au Au central 
collisions. A typical trigger particle {pt between 3.0 to 4.0 GeV/c) implies two associated regions in the considered 4> 
and r) ranges where the ridge particles lie. The ridge particles we consider are separated by 0.7 in r) {Ar)) so that one 
eliminates particles coming from jet production. The bulk of the ridge particles lie within 20° of the <p of the trigger 
particle. The width of the <f) spread is = 40°. 
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FIG. 10: On a trigger by trigger basis we accumulate particles from the ridge region (see Fig. 9). We then form a 
ratio the of pt spectrum for the ridge region particles divided by the pt spectrum for all particles in central Au Au 
events. We normalize the two pt spectra to have the same counts in the 2.0 GeV/c bin. This ratio or correlation 
for the bubbles which have a correlated pt among the partons (GFTM like) are the solid points. The random pt 
distribution among partons in the bubble are the open points. 
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y/SNN = 200 GeV requiring one trigger particle pt between 3 to 4 GeV/c and an associated particle pt above 2.0 
GeV/c. 



Triggered particles 




FIG. 12: The PBM generated CI correlation for the 0-5% centrality bin requiring one trigger particle pt above 3 
GeV/c and less than 4 GeV/c and another particle pt above 2.0 GeV/c plotted as a two dimensional A(j} vs. Arj 
perspective plot. The trigger requirements on this figure are the same as those on the experimental data in Fig. 11. 
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FIG. 13: The ridge signal is the piece of the CI correlation for the 0-5% centrality of Fig. 12 after removing all 
other particle pairs except the pairs that come from the same bubble. It is plotted as a two dimensional A</i vs. A77 
perspective plot. 



Jet Signal 




FIG. 14: The jet signal is left in the CI correlation after the contributions from the background and all the bubble 
particles are removed from the 0-5% centrality (with trigger requirements) of Fig. 12. It is plotted as a two dimensional 
A<^ vs. A77 perspective plot. 
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VII. STRONG CP VIOLATIONS OR 
CHERN-SIMONS TOPOLOGICAL CHARGE 

A. The Source FF. 

The strong CP problem remains one of the most 
outstanding puzzles of the Standard Model. Even 
though several possible solutions have been put for- 
ward it is not clear why CP invariance is respected 
by the strong interaction. It was shown however 
through a theorem by Vafa-Wittenfs?] that the true 
ground state of QCD cannot break CP. The part of 
the QCD Lagrangian that breaks CP is related to 
the gluon-gluon interaction term FF. The part of 
FF related to CP violations can be separated into 
a separate term which then can be varied by mul- 
tiplying this term by a parameter called 9. For the 
true QCD ground state 6 — (Vafa-Witten theo- 
rem). In the vicinity of the deconfined QCD vacuum 
metastable domains [ssj 9 non-zero could exist and 
not contradict the Vafa-Witten theorem since these 
are not the QCD ground state. The metastable do- 
mains CP phenomenon would manifest itself in spe- 
cific correlations of pion momenta ^3^ i3^]. 

B. Pionic measures of CP violation. 

The glasma flux tube model (GFTM)[l3] consid- 
ers the wavefunctions of the incoming projectiles, 
form sheets of CGCfTs'l at high energies that collide, 
interact, and evolve into high intensity color electric 
and magnetic fields. This collection of primordial 
fields is the Glasma[l^ [13]. Initially the Glasma 
is composed of only rapidity independent longitudi- 
nal (along the beam axis) color electric and mag- 
netic fields. These longitudinal color electric and 
magnetic fields generate topological Chern-Simons 
charge [21} through the FF term and becomes a 
source of CP violation. How much of these longi- 
tudinal color electric and magnetic fields are still 
present in the surface flux tubes when they have 
been pushed by the blast wave is a speculation of 
this paper for measuring strong CP violation? The 
color electric fleld which points along the flux tube 
axis causes an up quark to be accelerated in one di- 
rection along the beam axis, while the anti-up quark 
is accelerated in the other direction. So when a pair 
of quarks and anti-quarks are formed they separate 
along the beam axis leading to a separated tt'^ tt^ 
pair along this axis. The color magnetic field which 
also points along the flux tube axis (which is paral- 
lel to the beam axis) causes an up quark to rotate 
around the flux tube axis in one direction, while the 
anti-up quark is rotated in the other direction. So 
when a pair of quarks and anti-quarks are formed 



they will pickup or lose transverse momentum after 
the bubble is boosted radial outwards due to radial 
flow of the blast wave. These changes in pt will be 
transmitted to the 7r+ tt~ pairs. 

It is important to note that the CP violating asym- 
metries in TT"*" and 7r~ momenta arise through the 
Witten- Wess-Zumino^^ IsBj term. The quarks and 
anti quarks which, later form the 7r+ and , di- 
rectly respond to the color electric and color mag- 
netic fields and receive their boosts at the quark and 
anti quark level before they are pions. These boosts 
are transmitted to the pions at hadronization. Thus 
no external magnetic field is required in the method- 
ology followed in this paper. The following references 
demonstrate this:^ y, IH, . 

To represent the color electric field effect we as- 
sume as the first step for the results being shown 
in this paper; that we generate bubbles which have 
an added boost of 100 MeV/c to the quarks in the 
longitudinal momentum which represents the color 
electric effect. The 7r+ and the tt^ which form a pair 
are boosted in opposite directions along the beam 
axis. In a given bubble all the boosts are the same 
but vary in direction from bubble to bubble. For 
the color magnetic field effect, we give 100 MeV/c 
boosts to the transverse momentum in an opposite 
way to the tt"*" and the vr^ which form a pair. For 
each pair on one side of a bubble one of the charged 
pions boost is increased and the other is decreased. 
While on the other side of the bubble the pion for 
which the boost was increased is now decreased and 
the pion which was decreased is now increased by the 
100 MeV/c. All pairs for a given bubble are treated 
in the same way however each bubble is random on 
the sign of the pion which is chosen to be boosted 
on a given side. This addition to our model is used 
in the simulations of the following subsections. 



C. Color electric field pionic measure. 

Above we saw that pairs of positive and nega- 
tive pions should show a charge separation along the 
beam axis due to a boost in longitudinal momen- 
tum caused by the color electric field. A measure of 
this separation should be a difference in the pseudo- 
rapidity (At?) of the opposite sign pairs. This Ary 
measure has a well defined sign since we defined this 
difference measuring from the tt" to the 7r+. In or- 
der to form a correlation we must pick two pairs for 
comparison. The pairs have to come from the same 
bubble (a final state of an expanded flux tube) since 
we have shown by investigating the events generated 
by the model that pairs originating from different 
bubbles will not show this correlation. Therefore we 
require the 7r+ and the vr^ differ by 20° or less in </>. 
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Let us call the first pair A771. The next pair (A?72 
having the same requirement) has to also lie inside 
the same bubble to show this correlation. Thus we 
require that there is only 10° between the average (p 
of the pairs. This implies that at most in cf) no two 
pions from a given bubble can differ by more than 
30°. In Fig. 15 we show two pairs which would fall 
into the above cuts. At^i and A772 are positive in Fig. 
15. However if we would interchange the 7r+ and 7r~ 
on either pair the value of its A77 would change sign. 
Finally the mean value shown on Fig. 15 is the mid- 
point between the 7r+ and the 7r~ where one really 
uses the vector sum of the 7r+ and the tt" which 
moves this point toward the harder pion. 

Considering the above cuts we defined a correla- 
tion function where we combine pairs each having a 
At;. Our variable is related to the sum of the abso- 
lute values of the individual Ayy's (|A77i|-|-|A772|). We 
assign a sign to this sum such that if the sign of the 
individual Ayy's are the same it is a plus sign, while if 
they are different it is a minus sign. For the flux tube 
the color electric field extends over a large pseudora- 
pidity range therefore let us consider the separation 
of pairs IA77I greater than 0.9. For the numerator 
of the correlation function we consider all combina- 
tions of unique pairs (sign (|A?7i| -I- |A?72|)) from a 
given central Au Au event divided by a mixed event 
denominator created from pairs in different events. 
We determine the rescale of the mixed event denom- 
inator by considering the number of pairs of pairs for 
the case |A(^| lying between 50° and 60° for events 
and mixed events so that the overall ratio of this 
sample numerator to denominator is 1. By picking 
50° < |A(/)| < 60° we make sure we are not choosing 
pairs from the same bubble. For a simpler notation 
let (sign (|Aryi| -I- |A?72|)) = A?7i -f A?72 which varies 
from -4 to -|-4 since we have an over all 77 acceptance 
-1 to +1 (for the STAR TPC detector for which we 
calculated). The value being near ± 4 can happen 
when one has a hard pion with pt of 4 GeV/c (upper 
cut) at 77 = 1 with a soft pion pt of 0.8 GeV/c (lower 
cut) at 77 = -1 combined with another pair; a hard 
pion with pt of 4 GeV/c at 77 = -1 with a soft pion 
with Pt of 0.8 GeV/c at r/ = 1. 

In Fig. 16 we show the correlation function of op- 
posite sign charge-particle-pairs paired and binned 
by the variable A771 -I- Ar/2 with a cut |A?7| greater 
than 0.9 between the vector sums of the two pairs. 
The events are generated by the PBM@ and are 
charged particles of 0.8 < pt < 4.0 GeV/c, and I77I 

< 1, from Au Au collisions at ^/SNN = 200 GeV. 
Since we select pairs of pairs which are near each 
other in (j) they all together pick up the bubble cor- 
relation and thus these pairs of pairs over all show 
about a 4% correlation. In the 1.0 < Ar7i + Ari2 

< 2.0 region the correlation is 0.5% larger than the 



-2.0 < A771 + A772 < -1.0 region. This means there 
are more pairs of pairs aligned in the same direc- 
tion compared to pairs of pairs not aligned. This 
alignment is what is predicted by the color electric 
field effect presented above. In fact if one has plus 
minus pairs all aligned in the same direction and 
spread across a pseudorapidity range, locally at any 
place in the pseudorapidity range one would observe 
an increase of unlike sign charge pairs compared 
to like sign charge pairs. At small A77 unlike sign 
charge pairs are much larger than like sign charge 
pairs in both the PBM and the data which agree. 
See Figs. 10, 11, and 14 of Ref.Q. Figure 11 of 
Ref . j6i] compares the total correlation for unlike-sign 
charge pairs and like-sign charge pairs in the preci- 
sion STAR central production experiment for Au Au 
central collisions (0-10% centrality) at ^/snn—'^OO 
GeV, in the transverse momentum range 0.8 < pt 
< 2.0 GeV/c[l2|. The unlike-sign charge pairs are 
clearly larger in the region near A0 — A77 = 0.0. 
The increased correlation of the unlike-sign pairs is 
0.8% ± 0.002%. Figure 10 of Ref.0 shows that the 
PBM fit to these data gives the same results. Fig- 
ure 14 of Ref.Q shows that the CD = unlike-sign 
charge pairs minus like-sign charge pairs is positive 
for the experimental analysis. Therefore the unlike- 
sign charge pairs are considerably larger than the 
like-sign charge pairs. In fact this effect is so large 
and the alignment is so great that when one adds 
the unlike and like sign charge pairs correlations to- 
gether there is still a dip at small A77 and A^ see 
Fig. 4 and Fig. 7 of the present paper. The statisti- 
cal significance of this dip in the two high precision 
experiments done independently from different data 
sets gathered 2 years apart [H, [13] is huge. It would 
require a fluctuation of 14ct to remove the dip. This 
dip is also not due to any systematic error since both 
of the just cited precision experiments carefully in- 
vestigated that possibility; and found no evidence 
to challenge the reality of this dip. This highly sig- 
nificant dip (^14ct) means that like-sign pairs are 
removed as one approaches the region A?7 — Acj) — 
0.0 at a rate much larger than regular jet fragmen- 
tation. Jet fragmentation does not have a dip in its 
CI correlation. Thus this is very strong evidence for 
the predicted effect of the color electric field. 



D. Color magnetic pionic mesisure. 

After considering the color electric effect we turn 
to the color magnetic effect which causes up quarks 
to rotate around the flux tube axis in one direction, 
while the anti-up quarks rotate in the other direc- 
tion. So when a pair of quarks and anti-quarks are 
formed they will pickup or lose transverse momen- 
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FIG. 15: Pairs of pairs selected for forming a correlation due to the boost in longitudinal momentum caused by the 
color electric field. The largest separation we allow in (f> for plus minus pair 1 is A0i = 20°. This assures the pair is 
in the same bubble. The same is true for pair 2 so that it will also be contained in this same bubble. The mid-point 
for pair 1 and 2 represents the vector sum of pair 1 and 2 which moves toward the harder particle when the momenta 
differ. These mid-points can not be separated by more than lO" in A(j) in order to keep all four particles inside the 
same bubble since the correlation function is almost entirely generated within the same bubble. The A77 measure is 
the angle between the vector sum 1 compared to the vector sum 2 along the beam axis (for this case \Ari\ = 1.1). 
The positive sign for Ar/i comes from the fact that one moves in a positive 77 direction from negative to positive. The 
same is true for Ari2. If we would interchange the charge of the particles of the pairs the sign would change. 



turn. These changes in pt will be transmitted to 
the tt"*" and 7r~ pairs. It was previously shown in 
Sec. VII A that the CP violating asymmetries in the 
7r+ and the tt^ momenta arise through the Witten- 
Wess-Zumino term. The quarks and anti quarks 
which, later form the 7r+ and the tt", directly re- 
spond to the color electric and color magnetic fields 
and receive their boosts at the quark and anti quark 
level before they hadronize into pions. These boosts 
are transmitted to the 7r+ and the -k~ . 

In order to observe these differential pt changes 
one must select pairs on one side of the bubble in cj) 
and compare to other pairs on the other side of the 



same bubble which would lie around 40° to 48° away 
in (f). We defined a pair as a plus particle and minus 
particle with an opening angle (9) of 16° or less. We 
are also interested in pairs that are directly on the 
other side of the bubble. We require they are near in 
pseudorapidity {Ar] < 0.2). The above requirements 
constrain the four charged particles comprising both 
pairs to be contained in the same bubble and be close 
to being directly on opposite sides of the bubble (a 
final state expanded flux tube). The difference in pt 
changes due to and predicted by the color magnetic 
effect should give the tt"*" on one side of the bubble 
an increased pt and a decreased pt on the other side, 
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FIG. 16: Correlation lunction of pairs formed to exhibit the effects of longitudinal momentum boosts by the color 
electric field as defined in the text. This correlation shows that there are more aligned pairs (correlation is larger 
by ^0.5% between 1.0 < Ar?! + Ar]2 < 2.0 compared to between -2.0 < Ar]i + Ar]2 < -1.0). Arji + A772 is equal 
to |Ar?i| + |Ar72|. As explained in the text this means there are more pairs of pairs aligned in the same direction 
compared to pairs not aligned as predicted by the color electric field. The signs and more detail are also explained 
in the text. 



while for the tt" it will be the other way around. 
This will lead to an anti-alignment between pairs. 
In Fig. 17 wc show two pairs which would fall into 
the above cuts. Both pairs are at the limit of the 
opening angle cut 9i and 62 equal 16°. The pt of the 
plus particle for pair number 1 is 1.14 GcV/c, while 
the minus particle is 1.39 GeV/c. Thus APu is equal 
to -0.25 GeV/c. The pt of the plus particle for pair 
number 2 is 1.31 GeV/c, the minus particle is 0.91 
GeV/c and APta is equal to 0.40 GeV/c. Finally 
the mean value shown on Fig. 17 is the mid-point 
between the tt"*" and the tt~ where one really uses 
the vector sum of the 7r+ and the tt" which moves 
this point toward the harder pion. 

Considering the above cuts we defined a correla- 
tion function where we combine pairs each having a 



APf. Our variable is related to the sum of the abso- 
lute values of the individual APj's (| APn | + 1 AP(2 1)- 
We assign a sign to this sum such that if the sign of 
the individual AP^'s are the same it is a plus sign, 
while if they are different it is a minus sign. For 
the flux tube the color magnetic field extends over 
a large pseudorapidity range where quarks and anti- 
quarks rotate around the flux tube axis, therefore we 
want to sample pairs at the different sides of the tube 
making a separation in (p {A(p) between 40° to 48°. 
We are interested in sampling the pairs on the other 
side so we require the separation in r] {At]) be 0.2 or 
less. For the numerator of the correlation function 
we consider all combinations of unique pairs (sign 
(|APti| -t- |APt2|)) from a given central Au Au event 
divided by a mixed event denominator created from 
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FIG. 17: Pairs of pairs selected for forming a correlation exhibiting the effect of changes in pt due to the color 
magnetic field. The largest opening angle 8 for plus minus pairs is 16° or less. This opening angle assures that each 
pair has a high probability that it arises from a quark anti-quark pair. In this figure we have picked two pairs at 
this limit {Oi = 16° and 6*2 = 16°). The mid-point for pair 1 and 2 represents the vector sum of pair 1 and 2 which 
moves toward the harder particle when the momenta differ. These mid-points are chosen to have 40° < |A0| < 48° 
in order for the pairs to be on opposite sides of the bubble. Since we are interested in pairs directly across the bubble 
we make the At; separation be no more than 0.2. The difference in pt for pair 1 is APti = -0.25 GeV/c, while the 
difference in pt for pair 2 is APt2 = 0.40 GeV/c. The minus sign for 1 follows from the fact that the plus particle has 
1.14 GeV/c and the minus particle has 1.39 GeV/c. The plus sign for 2 follows from the fact that the plus particle 
has 1.31 GeV/c and the minus particle has 0.91 GeV/c. 



pairs in different events. We determine the rescale 
of the mixed event denominator by considering the 
number of pairs of pairs for the case |A?7| lying be- 
tween 1.2 and 1.5 plus any value of |A0| for events 
and mixed events so that the overall ratio of this 
sample numerator to denominator is 1. By picking 
this A?7 bin for all |A0| we have around the same 
pair count as the signal cut with the correlation 
of the bubbles being washed out. For a simpler no- 
tation let (sign (lAPfil + |APt2|)) = APn -f APt2 
which we plot in the range from -4 to -1-4 since we 
have an over all range 0.8 to 4.0 GeV/c. Thus the 



maximum magnitude of APf's is 3.2 GeV/c which 
makes APti -I- APt2 have a range of ±6.4. However 
the larger values near these range limits occur very 
rarely. 

In Fig. 18 we show the correlation function of op- 
posite sign charged-particle-pairs paired and binned 
by the variable APti -I- APt2 with a cut jAT^j less 
than 0.2 between the vector sums of the two pairs, 
and with 40° < |A(/)| < 48°. The events are gen- 
erated by the PBM[6j] and are charged particles of 
0.8 < Pt < 4.0 GeV/c, with < 1, from Au An 
collisions at y^SNN = 200 GeV. Since we select pairs 



22 



1.014 - 



1.012 



.01 



Correlation greater 0.8 and less 4.0 GeV/c 
Au Au Central 200 GeV 
lAril < 0.2 and 40^< |A^| < 48^ 



+ 



•••• 



• m 



• •• 



1.00S - 



1.006 - 



.004 



1 .002 

o 



0.998 - 



0.996 - 



AP„+AP,, 

FIG. 18: Correlation function of pairs of pairs formed for exhibiting the effect of changes in pt due to the color 
magnetic field as defined in the text. This correlation shows that there are more anti-aligned pairs. The correlation 
is larger by ~1.2% for APti + APt2 = -4.0 compared to APn + APt2 = 4.0. APti + APt2 is equal to | APti| + | APt2| 
where the sign is also explained in the text. 



of pairs which are near each other in (j) (40° to 48°) 
they all together pick up the bubble correlation and 
thus these pairs of pairs over all show about a 0.4% 
correlation. In the -4.0 < APti -I- APt2 < -1-0 region 
the correlation increases from 0.4% to 1%, while in 
the 1.0 < APti + APt2 < 4.0 region the correla- 
tion decreases from 0.4% to -0.2%. This means the 
pairs of pairs are anti-aligned at a higher rate than 
aligned. This anti-alignment is what is predicted by 
the color magnetic field effect presented above. In 
fact the anti-alignment increases with pt as the ratio 
bubble particle to background increases. However 
these predicted color magnetic effects have not been 
searched for yet and therefore there is no experimen- 
tal evidence for them. Experimental investigations 
of STAR TPC detector charged particle correlations 
provide a method of obtaining experimental confir- 
mation of the anti-alignment caused by the color 



magnetic field. 



VIII. SUMMARY AND DISCUSSION 

In this article we summarize the assumptions 
made and the reasoning that led to the develop- 
ment and construction of our parton bubble model 
(PBM)Q, which successfully explained the charge- 
particle-pair correlations in the central (0-10% cen- 
trality) ,JIWn = 200 GeV Au Au data[l2. The 
PBM is also consistent with the Au Au central col- 
lision HBT results. This is presented and discussed 
in Sec. 4 of Ref.[6j and Sec. I and II of this pa- 
per. In Sec. Ill we discussed extending our model, 
which was a central collision model, to be able to 
treat the geomery of bubble production for 0-80% 
collision centralities (PBME^lSj) such as measured 
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and analyzed in RHIC data[T3|. 

The bubbles represent a significant substructure of 
gluonic hot spots formed on the surface of a dense 
opaque fireball at kinetic freezeout. The origins of 
these hot spots come from a direct connection be- 
tween our Parton Bubble Model (PEMM) and the 
Glasma Flux Tube Model (GFTM)[i3|)- In the 
GFTM a flux tube is formed right after the initial 
collision of the Au Au system. This flux tube ex- 
tends over many units of pseudorapidity (A77). A 
blast wave gives the tubes near the surface trans- 
verse flow in the same way it gave flow to the bubbles 
in the PBM. This means the transverse momentum 
(pt) distribution of the flux tube is directly trans- 
lated to the Pt spectrum. In the PBM this flux tube 
is approximated by a sum of partons which are dis- 
tributed over this same large Arj region. Initially 
the transverse space is filled with flux tubes of large 
longitudinal extent but small transverse size ^Qj^. 
The flux tubes that are near the surface of the flre- 
ball get the largest radial flow and are emitted from 
the surface. As in the parton bubble model these 
partons shower and the higher pt particles escape 
the surface and do not interact. Qg is around 1 
GeV/c thus the size of the flux tube is about 1/4 
fm initially. The flux tubes near the surface are ini- 
tially at a radius ~5 fm. The </> angle wedge of the 
flux tube ~ 1/20 radians or ~3°. In Sec. V (also see 
Sec. IV) we connect the GFTM to the PBM: by as- 
suming the bubbles are the final state of a fiux tube 
at kinetic freezeout, and discussing evidence for this 
connection which is further developed in Sec. VI and 
Sec. VII. 

With the connection of the PBM to the GFTM 
two new predictions become possible for our PBM. 
The first is related to the fact that the blast wave 
radial flow given to the flux tube depends on where 
the tube is initially in the transverse plane of the 
colliding Au Au system. The tube gets the same 
radial boost all along its longitudinal length. This 
means that there is correlated pt among the partons 
of the bubble. In this paper we consider predictions 
we can make in regard to interesting topics and com- 
parisons with relevant data which exist by utilizing 
the parton bubble model (PBM [6]), and its features 
related to the glasma flux tube model fGFTMlflT}). 

Topic 1: The ridge is treated in Sec. VI. 

Topic 2: Strong CP violation (Chern-Simons 
topological charge) is treated in Sec. VII. 

We show in Sec. VI that if we trigger on particles 
with 3 to 4 GeV/c and correlate this trigger parti- 
cle with an other charged particle of greater than 
1.1 GeV/c, the PBM can produce a phenomenon 
very similar to the ridge[17i, 23 30]. (See Figs. 6- 
12). We then selected charged particles inside the 
ridge and predicted the correlation that one should 



observe when compared to the average charged par- 
ticles of the central Au Au collisions at ^snn = 200 
GeV[ll[l|. 

In Sec. VI D (Comparison to data): Triggered 
experimental angular correlations showing the ridge 
were presented at Quark Matter 2006 [3l|. Figure 7 
shows the experimental A(j) vs. Ar] CI correlation for 
0-10% central Au Au collisions at y^ijvw = 200; re- 
quiring one trigger particle pt between 3 to 4 GeV/c 
and an associated particle pt above 2.0 GeV/c. The 
yield is corrected for the finite A77 pair acceptance. 

For the PBM generator, we then form a two 
charged particle correlation between one charged 
particle with a pt between 3.0 to 4.0 GeV/c and 
another charged particle whose pt is greater than 
2.0 GeV/c. These are the same trigger conditions 
as in Ref.[3l| which is shown in Fig. 11, that 
shows the corrected pair yield in the central data. 
Fig. 12 shows the correlation function generated by 
the PBM which does not depend on the number of 
events analyzed. The two figures were shown to be in 
reasonable agreement when compared as explained 
previously in Sec. VI D. In Fig. 13 we show the 
ridge signal predicted by the PBM for very similar 
data but with 0-5% centrality. Figure 14 shows the 
extraction of the jet signal. Explanations are given 
in the text. 

The second prediction is a development of a pre- 
dictive pionic measure of the strong CP Violation. 
The GFTM flux tubes are made up of longitudi- 
nal color electric and magnetic fields which gener- 
ate topological Chern-Simons charge|21| through the 
FF term that becomes a source of CP violation. The 
color electric field which points along the fiux tube 
axis causes an up quark to be accelerated in one di- 
rection along the beam axis, while the anti-up quark 
is accelerated in the other direction. So when a pair 
of quarks and anti-quarks are formed they separate 
along the beam axis leading to a separated 7r+ and 
TT^ pair along this axis. The color magnetic field 
which also points along the flux tube axis (which 
is parallel to the beam axis) causes an up quark to 
rotate around the flux tube axis in one direction, 
while the anti-up quark rotates in the other direc- 
tion. So when a pair of quarks and anti-quarks are 
formed they will pickup or lose transverse momen- 
tum. These changes in pt will be transmitted to the 
7r+ and tt~ pairs. 

The above pionic measures of strong CP viola- 
tion are used to form correlation functions based 
on four particles composed of two pairs which are 
opposite sign charge-particle-pairs that are paired 
and binned. These four particle correlations accu- 
mulate from bubble to bubble by particles that are 
pushed or pulled (by the color electric field) and 
rotated (by the color magnetic field) in a right or 
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left handed direction. The longitudinal color electric 
field predicts aligned pairs in a pseudorapidity or A77 
measure. The longitudinal color magnetic field pre- 
dicts anti-aligned pairs in a transverse momentum 
or APf measure. The observations of these corre- 
lations would be a strong confirmation of this the- 
ory. The much larger unlike-sign pairs than like-sign 
pairs in the PBM and the data; and the strong dip of 
the CI correlation at small A77 (see Fig. 4 and Sec. 
VII C for full details) shows very strong evidence 
supporting the color electric alignment prediction in 
the y^ijvw — 200 GeV central Au Au collision data 
analyses at RHICfGl. [T2j. This highly significant dip 
(~14it) means that like-sign pairs are removed as 
one approaches the region Arj — A(j) = 0.0. Thus 
this is very strong evidence for the predicted effect 
of the color electric field. The color magnetic anti- 
aligned pairs in the transverse momentum predic- 
tion, treated in Sec. VII D as of now has not been 
observed or looked for. However our predicted spe- 
cific four charged particle correlations can be used to 
search for experimental evidence for the color mag- 
netic fields. 

Our success in demonstrating strong experimental 
evidence for the expected color electric field effects 
from previously published data suggests that the 
unique detailed correlations we have presented for 
searching for evidence for the predicted color mag- 
netic field effects should be urgently investigated. 



If we are lucky and the predicted color magnetic ef- 
fects can be confirmed experimentally we would have 
strong evidence for the following: 

1) CP is violated in the strong interaction in 
isolated local space time regions where topological 
charge [2 Ij is generated. 

2) The glasma flux tube model (GFTM) which 
was evolved from the color glass condensate (CGC) 
would be found to be consistent with a very signifi- 
cant experimental check. 

3) The parton bubble model event generator 
(PBM) is clearly closely connected to the CFTM. 
The bubble substructure strongly supported by the 
PBM is likely due to the final state of the fiux tube 
at kinetic freezeout. 
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